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ABSTRACT 

Recovery of the hypothetical supernova (SN) delay-time distribution (DTD) - the SN 
rate versus time that would follow a brief burst of star formation - can shed light 
on SN progenitors and physics, as well as on the timescales of chemical enrichment. 
Previous attempts to reconstruct the DTD have been based either on comparison 
of mean SN rates versus redshift to cosmic star- formation history (SFH), or on the 
comparison of SN rates among galaxies with different mean ages. Here, we present an 
approach to recover the SN DTD that avoids the averaging and loss of information 
of other schemes. We compare the SFHs of individual galaxies to the numbers of SNc 
discovered by a survey in each galaxy (generally zero, sometimes one SN, rarely a 
few). We apply the method to a subsample of 3505 galaxies, hosting 82 type-la SNc 
(SNe la) and 119 core-collapse supernovae (CC SNe), from the Lick Observatory SN 
Search (LOSS), that have SFHs reconstructed from Sloan Digital Sky Survey (SDSS) 
spectra. We find a > 2a SN la DTD signal in our shortest-delay, "prompt" bin at 

< 420 Myr. We identify and study a systematic error, due to the limited aperture of 
the SDSS spectroscopic fibres, that causes some of the prompt signal to leak to the 
later bins of the DTD. Taking this systematic error into account, we demonstrate that 
a prompt SN la contribution is required by the data at the > 99% confidence level. 
We further find a 4cr indication of SNe la that are "delayed" by > 2.4 Gyr. Thus, the 
data support the existence of both prompt and delayed SNe la. We measure the time 
integral over the SN DTD. For CC SNc we find a total yield of 0.010 ± 0.002 SNe 
per M Q formed, in excellent agreement with expectations, if all stars more massive 
than 8 M Q lead to visible SN explosions. This argues against scenarios in which the 
minimum mass for core-collapse SNe is > 10 M Q , or in which a significant fraction 
of massive stars collapse without an accompanying explosion. For SNe la, the time- 
integrated yield is 0.0023 ± 0.0006 SNe per M Q formed, most of them with delays 

< 2.4 Gyr. Finally, we show the robust performance of the method on simulated 
samples, and demonstrate that its application to already-existing SN samples, such as 
the full LOSS sample, but with complete and unbiased SFH estimates for the survey 
galaxies, could provide an accurate and detailed measurement of the SN la DTD. 
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Supernovae (SNe) figure prominently in many fields, 
whether in their roles as calibratable candles for cosmol- 



2 D. Maoz et al. 



ogy, as the major sources of intermediate-mass and heavy 
elements, as heaters of the interstellar medium, as accelera- 
tors of cosmic rays, and more. Physically, they are separated 
mainly into core-collapse SNe (CC SNe), which occur when 
the iron core of a massive star collapses to form a neutron 
star or a black hole, and type-la SNe (SNe la), which ex- 
plode when a degenerate carbon-oxygen stellar core, prob- 
ably a white dwarf (WD), approaches (or, rarely, exceeds) 
the Chandrasekhar mass, igniting the carbon and triggering 
a thermonuclear runaway. The two paths most often hy- 
pothesised for this mass growth in SNe la are the single- 
degenerate (SD) scenario, whereby a WD in a semidetached 
binary accretes matter from a main-sequence or evolved nor- 
mal companion star ( Whelan & Iben 1973) , and the double- 
degenerate (DD) scenario, in which two WDs merge (Iben & 
Tutukov 1984; Webbink 1984). Additional, less conventional, 
paths have also been considered (e.g., Tout 2005; Maoz & 
Mannucci 2008; Raskin et al. 2009a; Rosswog et al. 2009). 

For CC SNe, while many questions remain regarding 
the progenitors and the physics of particular subtypes, mas- 
sive progenitor stars have by now been identified in pre- 
explosion images in a growing number of cases (see Smartt 
2009, for a recent summary). This contrasts with the situa- 
tion for SNe la, where only one or two very ambiguous pro- 
genitor identifications exist (Voss & Nelemans 2008; Roelofs 
et al. 2008; Gonzalez Hernandez et al. 2009; Kerzendorf et 
al. 2009). We thus do not quite know what is exploding in 
a SN la, an unsatisfactory situation given the ubiquity and 
importance of these events. 

Driven by these problems, a major objective of SN stud- 
ies has been the recovery of the SN delay-time distribution 
(DTD). The DTD is the SN rate as a function of time that 
would be observed following a <5-function burst of star for- 
mation. (In other contexts, the DTD would be called the de- 
lay function, the transfer function, or the Green's function.) 
Knowledge of the DTD would be useful for understanding 
the route along which cosmic metal enrichment and energy 
input by SNe proceed, but no less important, for obtaining 
clues about the SN progenitor systems. Different progenitor 
stars, binary systems, and binary-evolution scenarios predict 
different DTDs. 

The lifetime of a star with the minimum initial mass 
that is thought to lead to a CC SN explosion, > 8 M©, sets 
a time division between CC SNe and SNe la in the DTD, 
at ~ 40 Myr. The precise mass border between core col- 
lapse and WD formation also depends on metallicity. Fur- 
thermore, initially lower-mass stars in tight binaries can be- 
come "rejuvenated" by mass transfer and explode as CC SNe 
somewhat later than this. 

For SNe la, the situation is much less clear. In both 
of the currently popular progenitor scenarios, SD and DD, 
calculations of the DTD depend on a series of assump- 
tions regarding initial conditions (initial mass function, bi- 
narity fraction, mass-ratio distribution, separation distri- 
bution), and complex physics (mass loss, mass transfer, 
common-envelope evolution, accretion) that is sometimes 
computationally intractable except in the most rudimentary, 
parametrized forms (e.g., Yungelson & Livio 2000; Hurley et 
al. 2002; Han & Podsiadlowski 2004; Nelemans et al. 2005; 
Ruiter et al. 2009; Bogomazov & Tutukov 2009; Meng & 
Yang 2010; Mennekens et al. 2010). In principle, observa- 
tional estimates of the DTD could rule out particular theo- 



retical models. Given the theoretical uncertainties, it is prob- 
ably more realistic that the observations simply provide a 
ground truth that successful models will need to reproduce. 

Previous attempts to recover the DTD have used SN 
rates measured in surveys of galaxies at different redshifts 
(i.e., different cosmic times), compared to cosmic star- 
formation histories (SFHs). This has been attempted for 
field surveys (Gal- Yam & Maoz 2004; Strolger et al. 2004; 
Poznanski et al. 2007; Dahlen et al. 2008) and galaxy-cluster 
surveys (Maoz & Gal- Yam 2004; Maoz et al. 2010). An al- 
ternative approach has been to look at the SN rates per 
unit stellar mass in galaxies of particular types (star form- 
ing, quiescent, etc.), and to attempt to assign to each type 
a "formation age," or some generic, simple, SFH (e.g., Man- 
nucci et al. 2005, 2006; Sullivan et al. 2006; Totani et al. 
2008; Pritchet et al. 2008; Raskin et al. 2009b). 

Results have been controversial and often contradictory. 
For example, Dahlen et al. (2004, 2008) have argued for a 
SN la DTD that is peaked at a delay of ~ 3 Gyr, with few 
SNe la at delays that are much shorter or longer. In con- 
trast, Mannucci et al. (2005, 2006), Scannapieco & Bildsten 
(2005), and Sullivan et al. (2006) have found evidence for 
the existence of comparable numbers of both "prompt" and 
"delayed" SNe la: the former explode within ~ 500 Myr 
(or perhaps even within 100 Myr) of star formatiorQ, while 
the latter may have delays as long as 10 Gyr. The SN la 
rate has been described as the sum of two components, 
one proportional to stellar mass and the other proportional 
to the CC SN rate (Mannucci et al. 2005). In the simi- 
lar U A + B" parameterisation introduced by Scannapieco 
& Bildsten (2005), the prompt-component rate is propor- 
tional to the star- formation rate (SFR). The two SN la 
components need not represent two distinct physical pop- 
ulations. Instead, they could constitute the SNe included in 
two coarsely sampled time bins of what is in reality a contin- 
uous DTD. For example, Pritchet et al. (2008) have argued 
that a £-°' 5±0 - 2 power-law DTD provides an improved fit, 
compared to the A + B model, to the dependence of SN 
rates on galaxy SFR, as measured in the Supernova Legacy 
Survey. We also note that a truly bimodal DTD, if it exists, 
could arise either from two different coexisting progenitor 
paths (e.g., DD and SD), or from bimodality in some sec- 
ondary parameter, such as the binary separation distribu- 
tion, of a single explosion mechanism. Regardless, there is 
currently an unclear picture on the form of the DTD, even 
at the most coarse resolution level. 

A shortcoming of the approaches described above for 
recovering the DTD is that they involve averaging over the 
galaxy population (i.e., all the SNe are assumed to come 
from the entire host population considered), or averaging 
over time (i.e., the detailed SFH of a galaxy is represented 
by a single "age" or simplified history for all galaxies of a 
certain type). Consequently, these approaches involve loss 

1 We note that diverse delay ranges have been associated in the 
literature with the term "prompt" — e.g., < 100 Myr (Mannucci 
et al. 2006); < 180 Myr (Aubourg et al. 2008); 200-500 Myr 
(Raskin et al. 2009b); < 350, < 700 Myr, or < 1 Gyr (Scanna- 
pieco & Bildsten 2005). The term therefore generally labels delays 
of roughly a few hundred Myr. In our analysis of the Lick Obser- 
vatory SN Search herein, we will define prompt SNe la as those 
with delays < 420 Myr. 
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of information, and potential systematic errors (e.g., due to 
unrepresentative simplified histories). 

In this paper, the fourth in a series analysing SN rates 
from the Lick Observatory SN Search (LOSS), we introduce 
a new approach to recover the DTD, by posing DTD in- 
version as a discretised linear problem. In this process, we 
use all of the available information on the SFHs of individ- 
ual galaxies, rather than averaging rates over many types 
of galaxies in some redshift interval, or assigning a mean 
star-formation-weighted age to each galaxy. In §2 below, 
we present the method, and in §3, we apply the method 
to a subsample of LOSS galaxies and their SNe. We test the 
method's performance on simulated SN surveys in §4. Our 
results are summarised and we discuss some future prospects 
in §5. 

2 RECONSTRUCTION OF THE SN DTD — 
METHOD 

Consider a sample of N galaxies. The SN rate in the ith 
galaxy observed at cosmic time t is given by the convolution 

n(t) = / Sj(t-T)*(T)dT, (l) 
Jo 

where Si(t) is the SFR versus cosmic time of the ith galaxy 
(stellar mass formed per unit time), *P(r) is the DTD (SNe 
per unit time per unit stellar mass formed), and the integra- 
tion is from the Big Bang (t = 0) to the time of observation. 
For the purpose of this paper, we assume that the DTD is a 
universal function: it is the same in all galaxies, independent 
of environment, metallicity, and cosmic time — a simplifying 
assumption that may be invalid at some level. For example, 
a dependence of SN delay time on metallicity is expected in 
some models (e.g., Kobayashi et al. 2000). Similarly, varia- 
tions in the initial mass function (IMF) with cosmic times 
or environment would also lead to a variable DTD, but we 
will again ignore this possibility in the present context. 

In contrast to the averaging approaches followed in the 
past (see §1), we will attempt to recover the DTD by directly 
inverting a linear, discretized version of Eq. [1] where the de- 
tailed history of every individual galaxy or galaxy subunit is 
taken into account. Suppose the SFHs of the i — 1,2, N 
galaxies monitored as part of a SN survey are known (e.g., 
based on reconstruction of their stellar populations), with 
a temporal resolution that permits binning the stellar mass 
formed in each galaxy into j = 1, 2, .., K discrete time bins, 
where increasing j corresponds to increasing lookback time. 
The time bins need not necessarily be equal, and generally 
will not be, since the temporal resolution of the SFH recon- 
struction degrades with increasing lookback time. For the 
ith galaxy in the survey, the stellar mass formed in the jih 
time bin is rriij. The mean of the DTD over the jih bin 
(corresponding to a delay range equal to the lookback-time 
range of the jih bin in the SFH) is tyj. Then the integration 
in Eq. [T]can be approximated as a sum, 

K 

r. y\/ ( . ; >l<,. (2) 

3 = 1 

where n, the SN rate in a given galaxy, is measured at a 
particular cosmic time (e.g., corresponding to the redshift 
of the particular SN survey). Given a survey of N galaxies, 



each with an observed SN rate, n, and a known binned 
SFH, rriij, one could, in principle, algebraically invert this 
set of linear equations and recover the best-fit parameters 
describing the binned DTD: * = (*i, * 2 , — , ^k)- 

In practice, on human timescales SNe in a given galaxy 
are rare events (n -C 1 yr^ 1 ). Supernova surveys therefore 
monitor many galaxies, and record the number of SNe dis- 
covered in every galaxy. For a given model DTD, 4', the ith 
galaxy will have an expected number of SNe 

A* = nU, (3) 

where ti is the effective visibility time (often called the "con- 
trol time") during which a SN of a particular type would 
have been visible (given the actual on-target monitoring 
time, the distance to the galaxy, the flux limits of the survey, 
and the detection efficiency). Since Xi <C 1, the number of 
SNe observed in the ith galaxy, m, obeys a Poisson proba- 
bility distribution with expectation value Xi, 

P(ni\Xi) = (e-M^M!, (4) 

where m is for most of the galaxies, 1 for some of the 
galaxies, and more than 1 for very few galaxies. 

2.1 Maximum-likelihood DTD recovery 

We now introduce a nonparametric, maximum-likelihood 
method to recover the DTD and its uncertainties. Consid- 
ering a set of model DTDs, the likelihood of a particular 
DTD, given the set of measurements m, tin, is 

N 

L = ]JP(n l \X l ). (5) 

i=l 

More conveniently, the log of the likelihood is 

N N N 

InL = J2lnP(n i \X i ) = - ^A« + ^ln(A? i /n i !), (6) 

z— 1 z — 1 i — 1 

where obviously only galaxies hosting SNe contribute to the 
second term. The best-fitting model can be found by scan- 
ning the parameter space of the vector \l/ for the value that 
maximizes the log-likelihood. This procedure naturally al- 
lows restricting the DTD to have only positive values, as 
physically required (a negative SN rate is meaningless). 

The covariance matrix Cjk of the uncertainties in the 
best-fit parameters can be found (e.g., Press et al. f992) by 
calculating the curvature matrix, 

_ 1 d 2 In L _ y> d[lnP(ni\Xi)] d{ln P(m\Xi)} 

J I— 1 J 

N 

= ^2t 2 (ni/Xi - l) 2 mijm ik , (7) 

!=1 

and inverting it, 

[C] = [a]' 1 . (8) 

Because the values of the DTD are constrained to be posi- 
tive, if the maximum-likelihood value of a DTD component, 
tyj, is close to zero, the square root of its variance, \fCjj, 
will not represent well its la uncertainty range. An alterna- 
tive, more reliable, procedure is to perform a Monte-Carlo 
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simulation in which many mock surveys are produced, each 
having the same galaxies, SFHs, and visibility times as the 
real survey, and having expectation values Xi based on the 
best-fit DTD, but with the number of SNe in every galaxy, 
m, drawn from a Poisson distribution according to A;. The 
maximum-likelihood DTD, vE", is found for every realization. 
From the distribution of the values of every component, vp^, 
over all the realizations, one can estimate the range encom- 
passing, say, ±34% of the cases. 

The above approach for recovering the DTD has several 
advantages over previous methods. First, all of the known 
information in the survey is included in the analysis in a sta- 
tistically rigourous way, including the fact that many (usu- 
ally most) of the galaxies did not host any SNe. Further- 
more, the calculation is easily generalized to cases where 
the galaxies are not all at the same distances (e.g., combi- 
nations of surveys done at different redshifts) - one simply 
needs to use the appropriate SFH bins and visibility times 
for every galaxy. In fact, assuming that the DTD is a uni- 
versal function, it is straightforward to include in a single 
analysis the data from completely disparate SN surveys. For 
example, one could combine the results of normal SN sur- 
veys with unconventional SN "surveys," in which the SN 
rate is measured based on SN remnants in small subunits of 
a few nearby galaxies (Maoz & Badenes 2010). 

The number and resolution of the time bins used in the 
analysis will naturally depend on the quality of the data. 
Larger numbers of observed SNe, JVtot , as well as better data 
on the parent stellar populations (integrated colors and/or 
spectra), will permit a larger number of independent SFH 
and DTD time bins, and will thus improve the time resolu- 
tion of the recovered DTD. We quantify this in §4, below. 



3 APPLICATION TO THE LOSS-SDSS 
SAMPLE 

We now apply our method to the SN survey data obtained 
by considering all LOSS galaxies with SFHs based on spec- 
troscopy from the Sloan Digital Sky Survey (SDSS). First, 
we summarize briefly the essentials of each of these surveys 
and of the sample resulting from their intersection. 

3.1 The Lick Observatory SN Search 

The LOSS is an ongoing survey for SNe in a sample of 
~15,000 nearby (redshift z < 0.05) galaxies, conducted with 
the Katzman Automatic Imaging Telescope (KAIT) at Lick 
Observatory (Li et al. 2000; Filippenko et al. 2001; Filip- 
penko et al. 2010). KAIT is a fully robotic telescope whose 
control system checks the weather and performs observa- 
tions with a dedicated CCD camera without human inter- 
vention. The data are automatically processed through an 
image-subtraction pipeline, and candidate SNe are flagged 
and visually inspected. The promising SN candidates are re- 
observed and the confirmed SNe are reported to the Central 
Bureau of Astronomical Telegrams. 

A series of papers (Filippenko et al. 2010; Leaman et al. 
2010; Li et al. 2010a,b) present the details and first results 
on the rates of SNe in the local Universe based on LOSS, 
using a sample of 1036 SNe discovered in more than 2 mil- 
lion observations between March 1998 through the end of 



2008. This is the largest and most homogeneous set of SN 
statistics ever assembled for the determination of local SN 
rates. Filippenko et al. (2010) describe the instrumentation 
and the details of the SN survey. Leaman et al. (2010, Paper 
I in this series) present the control-time calculation for the 
galaxies in the sample, for SNe of different types and their 
luminosity functions (LFs), and details of the galaxy and 
SN samples used in the rate calculations. Monte-Carlo sim- 
ulations are used to determine the limiting magnitude and 
the SN detection efficiency in each LOSS search image. 

Li et al. (2010a, Paper II in the series) discuss the ob- 
served SN LF using a volume-limited sample (D < 60 Mpc 
for CC SNe and D < 80 Mpc for SNe la) of 177 SNe, each 
with detailed spectroscopic classification and peak magni- 
tude from dedicated photometric follow-up images or unfil- 
tered survey images. These observed LFs solve two issues 
that have plagued historical SN-rate calculations - the in- 
trinsic luminosity distribution of SNe and the host-galaxy 
extinction. Finally, Li et al. (2010b, Paper III) combine all 
of the above ingredients, obtaining control times for different 
types of SNe and for each galaxy, based on its monitoring 
history, the observed LFs, and the limiting magnitudes and 
detection efficiencies of the search images. These are used to 
derive SN rates for SNe of different types, as a function of 
various galaxy properties. 

3.2 VESPA star-formation histories of SDSS 
galaxies 

The SDSS (York et al. 2000) is a survey of ~ 10 4 deg 2 of the 
north Galactic cap, consisting of imaging in five photometric 
bands (u, g,r,i, z), and 3"-aperture fibre spectroscopy of ~ 
10 6 targets, mostly galaxies, with r < 18 mag. Tojeiro et al. 
(2009) performed spectral synthesis modeling of all galaxy 
spectra in the SDSS using their VESPA code (Tojeiro et al. 
2007). VESPA uses all of the available absorption features, 
as well as the shape of the continuum, to deconvolve the 
observed spectra and obtain an estimate of the SFH. In order 
to recover the maximum amount of reliable information, the 
number of time bins used is variable, and depends on the 
quality of the data on each galaxy. At the highest resolution, 
VESPA uses 16 age bins, logarithmically spaced between 
0.002 Gyr and to = 13.7 Gyr, the age of the Universe. When 
data do not have sufficiently high signal-to-noise ratio for a 
fully resolved reconstruction, pairs of adjacent time bins are 
averaged. This process may be repeated down to the last 
two remaining bins. In the end, the SFH of every galaxy is 
computed using a different set of time bins. 

VESPA masses are calculated assuming a Kroupa 
(2007) IMF; all of our results will therefore include this as- 
sumption implicitly. Bell et al. (2003) have shown that the 
Kroupa IMF gives a similar total stellar mass to that of a 
"diet Salpeter" IMF, obtained by multiplying by 0.7 the to- 
tal mass of the original Salpeter (1955) IMF, to account for 
the reduced number of low-mass stars. Thus, our results will 
be comparable to other SN rate studies, such as Mannucci et 
al. (2005, 2006), that have assumed the diet-Salpeter IMF. 

3.3 The LOSS-SDSS- VESPA sample 

With the kind assistance of R. Tojeiro, we have derived 
VESPA SFHs for all of the LOSS galaxies that have SDSS 
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spectroscopy. Our main sample of LOSS galaxies with 
VESPA SFHs consists of 3505 galaxies that hosted 201 SNe, 
among them 82 SNe la, 93 SNe II, and 26 SNe Ibc (see 
Filippenko 1997 for a review of SN types; here we classify 
SNe lb and Ic as "Ibc"). An alternative sample, using a dif- 
ferent VESPA dust model (see below), has slightly different 
numbers: 3508 galaxies hosting 202 SNe, among them 81 
SNe la, 94 SNe II, and 27 SNe Ibc. These numbers of SNe 
are consistent with those expected based on the fraction of 
the total LOSS visibility time that is in the VESPA galaxies. 
Other than the increase in Poisson errors due to the smaller 
number of SNe, the limitation of the calculation to VESPA 
galaxies should have no effect on our DTD reconstruction. 

For the majority of these galaxies, it is practical to sep- 
arate the SFHs into no more than four time bins: — 70 Myr, 
70 - 420 Myr, 420 Myr-2.4 Gyr, and > 2.4 Gyr. These cor- 
respond to the bins labeled 24, 25, 26, and 27 in Tojeiro et 
al. (2009) . However, for each galaxy, the distribution of mass 
between the first two bins depends strongly on the partic- 
ular dust model that is assumed, as well as on the spectral 
synthesis code that is used - Bruzual & Chariot (2003), or 
Maraston (2005) (R. Tojeiro, 2009, private communication). 
In fact, for some models, only a small fraction of the galaxies 
have any mass formed in the second time bin, a problem that 
persists at some level in other models as well. Due to this de- 
generacy, we have chosen to combine the first two time bins 
into a single bin, of — 420 Myr. As the choice of stellar pop- 
ulation model parameters affects the results in other time 
bins as well, we have used two alternative VESPA SFH re- 
constructions. One assumes a single dust component for each 
galaxy, while the other allows separate dust components for 
the young and the old stellar populations, effectively intro- 
ducing an additional free parameter to the modeling (see 
Tojeiro et al. 2009, for details). Both reconstructions use 
the Maraston (2005) spectral synthesis models. Using these 
different SFH reconstructions gives some idea of the sys- 
tematic errors in the DTD reconstruction arising from the 
uncertainty in the SFH (until now, we have treated the SFH 
values, rriij, as error-free independent variables). 

3.4 The core-collapse-SN DTD 

We begin by deriving the DTD of the CC SNe in the LOSS- 
SDSS sample. Core-collapse SNe explode within < 40 Myr 
of star formation, and therefore their DTD should have zero 
amplitude on timescales much longer than this. This can 
provide a first test of our sample and of the DTD recovery 
method. 



3.4-1 Full-sample core-collapse SN DTD 

Figure [1] shows the reconstructed CC SN DTD, and its 
uncertainties from the Monte-Carlo simulations described 
above, for the full sample of 3505 LOSS galaxies with 
VESPA SFHs using a single dust parameter. The SNe in 
this sample consist of 93 SNe II and 26 SNe lb and Ic. 
Since SNe II and SNe Ibc have different visibility times, we 
have derived the DTD for each subsample separately, added 
the resulting DTD amplitudes in the corresponding bins, 
and added the uncertainties in quadrature. Like the SFHs, 
the DTD in Figure [T] has three time bins, corresponding to 



0.1 


: 






119 CC SNe in 


3505 


galaxies — 








• 


1 — parameter 


dust model : 




I 


true CC 


signal 






















30.01 








1 #— < 






"0 








leaked signal 






t 0.001 
u 

!>, 














Z 
in 














0.0001 












i 


10- 6 
























10 10 






10 8 




10 9 





Delay time [yr] 

Figure 1. Best-fit delay-time distribution found for the full 
LOSS-SDSS core-collapse SN sample using the VESPA SFH re- 
construction with a single dust parameter. Points mark the best- 
fit values for each time bin, whose range is indicated by the hor- 
izontal error bars. Vertical error bars show the most likely 68% 
range. Note the nonzero amplitude of the DTD in the delayed 
bins, due to "leakage" from the prompt bin resulting from incor- 
rect characterization of the full stellar populations of the galaxies 
by the limited aperture of the SDSS spectra. 
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Figure 2. Same as Fig. [T] but using a subsample obtained by 
culling the LOSS-SDSS sample of all galaxies of Hubble type Sa- 
Sbc, for which SDSS fibre spectra do not represent well the the 
full stellar population of a galaxy. Symbols are as in Fig. [T] Note 
the enhanced ratio of rates between the first and second bins, 
due to reduced leakage. The best-fit rate in the third bin is zero, 
and its 95% confidence upper limit is marked. Nevertheless, the 
leakage is not eliminated completely, and the statistical errors are 
larger because of the reduced sample. 

< 420 Myr ("prompt" SNe), 420 Myr - 2.4 Gyr ("medium 
delay" SNe), and > 2.4 Gyr ("delayed" SNe). The best-fit 
values and the uncertainties for these DTDs, as well as ad- 
ditional ones discussed below, are listed in Table [1] 

The prompt, < 420 Myr bin, in which all CC SNe should 
reside, indeed has a clear, > 7a signa0. However, there is 

2 We will henceforth use the phrasing u Xa signal" to denote the 
number of —la (—34%) errors by which a best-fit DTD level is 
above zero. Naturally, in some cases this will not correspond to 
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Table 1. Summary of DTD reconstructions 



Sample 




N S N 




*2 


*3 


N SN /M 








0-0.42 Gyr 


0.42-2.4 Gyr 


2.4-14 Gyr 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


CC SNe 


Full CC 


3505 


119 


770 ± 100 


127 ± 35 


8.6 ± 8.5 


0.00585 ± 0.00085 


No Sa-Sbc 


1900 


50 


1320 ± 230 


101 ± 66 


< 5.3 


0.00750 ± 0.00150 


D < 100 Mpc 


1951 


92 


1037 ± 146 


192 ± 62 


< 4.3 


0.00920 ± 0.00130 


D < 65 Mpc 


851 


45 


1430 ± 340 


223 ± 165 


< 13.4 


0.01040 ± 0.00360 



SNe la 



Full, 1-dust 3505 82 84±|| 37+^ 2.6±";| 0.00137 ± 0.00032 

Full, 2-dust 3508 81 26±2o 20 -i2 3 - 3 -o.t 0.00086 ± 0.00024 

No Sa-Sbc 1900 49 136+^° 55^ ' 3 -°-o.l 0.00200 ± 0.00060 

Column header explanations: 

(1) - Sample used to derive DTD: Full CC - full core-collapse (types II and Ibc) SN sample with VESPA SFH 
reconstructions using a single-parameter dust model; No Sa-Sbc — sample excluding all galaxies of Hubble type 
Sa through Sbc, to avoid small-fibre aperture effects; D < 100 Mpc - core-collapse SN sample in galaxies within 
100 Mpc: D < 65 Mpc - core-collapse SN sample in galaxies within 65 Mpc. Full, 1-dust - full SN la sample, with 
VESPA SFH reconstructions using a single-parameter dust model; Full, 2-dust — full SN la sample, with VESPA 
SFH reconstructions using a two-parameter dust model. 

(2) - Number of galaxies in sample. 

(3) - Number of SNe in sample. 

(4-6) - DTD rates and 68% uncertainty ranges, in units of 10 — 4 SNe yr — 1 (10 10 Mq)" 1 . Upper limits, where given, 
correspond to a best fit of 0, and to the 95% confidence limit. 
(7) - Time-integrated DTD, in units of SNe M" 1 . 



also a strong 3.5cr-level DTD signal in the 420 Myr-2.4 Gyr 
bin, where no CC explosions are expected. As we will demon- 
strate and quantify in i]3.4.2l below. this result is due to the 
following systematic error. The SDSS spectra of many of the 
LOSS galaxies (which are relatively nearby, and hence large 
in angle) are dominated by the old populations in the cen- 
tre of each galaxy, due to the limited 3" aperture size of the 
SDSS fibres. However, many of the CC SNe explode in the 
outer regions, where there is ongoing star formation that is 
invisible to the SDSS spectroscopy. As a result, our DTD 
solution mistakenly associates these CC SNe with an old 
population, and hence a large delay. We also note that the 
integrals over these two bins of the best-fit DTD are com- 
parable: *iAii = 0.00325 and * 2 At 2 = 0.00251 CC SNe 
per formed, respectively, where Ati = 420 Myr and 
At 2 = 1-98 Gyr are the lengths of the time bins. Thus, 
about 40% of the prompt-bin signal has leaked into bin 2. It 
is unavoidable that the same process will distort the results 
for the DTD of SNe la in the same way. 

In the third bin, > 2.4 Gyr, the best-fit amplitude is 
again nonzero, but only at the la level, and the integral 
over this bin, ^^Ats = 9 x 10 -5 Mq , is just a few percent 
of the first two bins. Thus, there is little leakage of power to 
bin 3. 



the Gaussian probability associated with Xcr, either for obtaining 
zero given the best-fit result, or for obtaining the best-fit result 
given zero. We will discuss these subtleties where relevant. 



3.4-2 Subsamples selected to reduce crosstalk 

To understand the effect of the small-aperture SDSS fibres, 
and in an attempt to find a method of obtaining cleaner 
signals in the DTD, we have experimented with culling the 
galaxy sample in various ways, based on the galaxy proper- 
ties. For example, we expect that much of the "leak" that 
we see is caused by early-type spirals, in which the spectrum 
of the bulge, covered by the SDSS fibre, is highly unrepre- 
sentative of the galaxy disk. In contrast, elliptical galaxies 
and late-type spirals both have more spatially uniform stel- 
lar populations. We have therefore derived the DTD, shown 
in Figure [2] for CC SNe from a subsample that excludes all 
LOSS galaxies with Hubble types 3 to 5 in the notation of 
Leaman et al. (2010), corresponding approximately to types 
Sa to Sbc. This reduces the sample size to 1900 galaxies, 
hosting 34 SNe II and 16 SNe Ibc (see Table 0. As ex- 
pected, the prompt signal in is strengthened by a factor 
2, while ^2 is reduced in amplitude. The best-fit value for 
^3 is reduced to 0. 

Unfortunately, as seen in Figure [2] and in the numbers 
in Table[T] hand-in-hand with the increase/decrease in DTD 
amplitude of the first/second bin, the smaller SN numbers 
in the reduced sample lead to larger statistical errors. We 
have experimented with additional sample culling methods. 
For example, the SDSS database gives, for each galaxy, the 
magnitude in each photometric band in an aperture of the 
size of the spectrograph fibre, and various measures of the 
total magnitude of the galaxy. One can thus select to in- 
clude in the DTD analysis only galaxies with no or weak 
radial color gradients. Additional or alternative cuts can be 
made based on galaxy size, distance, or fraction of total light 
within the fibre aperture (i.e., degree of concentration). We 
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LOSS-SDSS best-tit DDT 
82 SNe la in 3505 galaxies 
1— parameter dust model 
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Figure 3. Best-fit SN la delay-time distribution found for the 
LOSS-SDSS sample using the VESPA SFH reconstruction with 
a single dust parameter. Points mark the best-fit values for each 
time bin, whose range is indicated by the horizontal error bars. 
Solid vertical error bars show the most-likely 68% range, based 
on Monte-Carlo simulations that use the best-fit DTD values. 
Dashed vertical error bars (slightly shifted to the right, for clarity) 
show the Gaussian la errors from calculation of the covariance 
of the parameters. The Gaussian errors are symmetric about the 
best fit, but appear asymmetric in the plot because of the loga- 
rithmic scale. 



LOSS-SDSS best-fit DDT 
81 SNe la in 350B galaxies 
2— parameter dust model 
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Figure 4. Same as Fig. [3] but using a VESPA SFH reconstruction 
permitting separate dust components for the young and the old 
stellar populations of each galaxy. 



find that, as with the case of the selection (above) according 
to Hubble type, the more stringent the selection, the lower 
the leak out of the first bin in the CC SN DTD. At the same 
time, the lowered sample size increases the statistical errors. 

3.5 The SN la DTD 

Keeping in mind the systematics evidenced in the case of 
the CC SNe, we now proceed to derive the DTD for the 
SNe la in the sample. Figures [3HS] show the reconstructed 
SN la DTDs, and their uncertainties, for the LOSS-SDSS 
sample, using the two different VESPA-based SFHs. As for 
the CC SNe, the best-fit values and the uncertainties for 
these DTDs are listed in Table □ 



3.5.1 Full-sample SN la DTD, one-parameter dust model 

In the DTD based on the SFHs with a one-parameter dust 
model, we obtain, in the "prompt" (< 420 Myr) bin, a best- 
fit value of *i = 0.0084±°;°^ 4 SNe yr -1 (10 10 M ) -1 . (In 
other words, following a 5-function starburst that forms a 
stellar mass of 10 10 M©, the mean SN la rate in this stellar 
population over the first 420 Myr will be 0.0084±g;^|| SNe 
yr -1 .) This result implies a > 2a detection of a prompt SN la 
component. As before, the uncertainties we quote are the 
statistical errors, as derived from the most probable ±34% 
range in Monte-Carlo simulations, using the best-fit values 
as input. A model with the above best-fit value of as 
input yields a recovered of in fewer than 3% of the 
Monte-Carlo realizations. Conversely, a model with an input 
value of ^fi = yields a recovered ^i ^ 0.0084 (in the same 
units as before) in < 0.2% of the realizations. Models with 
an input value of *i < 0.0031 yield a recovered *i > 0.0084 
in < 5% of Monte-Carlo realizations. Thus = 0.0031 can 
be considered a 95%-confidence lower limit on the level of a 
prompt SN la component. 

As another test of the presence of a prompt SN la com- 
ponent, we have forced = in the DTD reconstruction. 
Compared to the previous result, the log of the likelihood 
of the best-fit model decreases by 3.0. Since 2AlnL « A\ 2 , 
this indicates that a prompt component in the model im- 
proves the fit to the data at the 2.5a level. Thus, this test 
supports, at the ~ 99% confidence level, the existence of 
SNe la that explode within 420 Myr after star formation. 

Since the uncertainties here are larger than was the 
case with the CC SN DTD, we can use this case to com- 
pare the adequacy of the errors calculated directly using the 
covariance-matrix formalism to those from the Monte-Carlo 
simulations. From Figs. [3HU we see that the analytic la 
errors generally match well the errors estimated from the 
simulations. From these simulations, we also find that the 
likelihood of the best-fit model is attained in at least 10% 
of simulated trials, indicating that the best-fit model is ac- 
ceptable in an absolute sense as well. 

The mean rate found above in the first bin, ^i, can 
be translated into an equivalent value of the "B" parame- 
ter (Scannapieco & Bildsten 2005), the constant of propor- 
tionality between the SFR and the "prompt" SN la rate 
(i.e., B is the number of prompt SNe per unit stellar mass 
formed). Multiplying *Pi by 420 Myr to integrate over the 
bin, and dividing by 10 10 to express the result per solar 
mass, the best-fit DTD value in the prompt bin implies 
B « (3.5 ± 1.7) x 10 _4 Mq 1 (or B < 8 x lO^M" 1 for 
the 2a upper limit). This is several times lower than the 
values of B estimated by studies that compare the SN la 
rate per unit mass and the SFR in blue, vigorously star- 
forming galaxies [see recent summary and intercomparison 
in Maoz (2008); observables that are quoted there per unit 
stellar mass formed need to be divided by 0.7 to convert 
from the pure Salpeter IMF assumption to the low-mass- 
truncated or "diet" IMFs considered here - Kroupa (2007) 
or, equivalently in terms of integrated mass, diet-Salpeter 
(Bell et al. 2003); see §3.2]. Such studies have usually found 
values of B = (1 — 3) x 10~ 3 Mg 1 (an exception is the value 
B = (3.9 ± 0.7) x 10 -4 M" 1 found by Sullivan et a. 2006). 

As we showed in ^3.4.21 above, at least part of this dis- 
crepancy must be due to the limited 3" aperture size of the 
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SDSS fibres. Just as our CC SN DTD solution mistakenly 
associated CC SNe with an old population, and hence a large 
delay, it will do so for prompt SNe la that are associated with 
star-forming regions in the outer parts of a galaxy, outside 
the bulge-light-dominated fibre aperture. In fact, the leak 
should be comparable to the 40% effect we saw for CC SNe. 
Thus, the true 'I'i level in Figure [3] is likely about double 
the derived, "leaky" one, and thus coincident, within the 
uncertainties, with the B-parameter values found by other 
studies, above. By the same token, our > 2a significance on 
the detection of a prompt component is a lower limit. 

Another contribution to the difference between the 
time-integrated value and published B values may be 
that the B parameter found from SN rate measurements in 
galaxies with the highest SFRs are tracing a SN la popula- 
tion with an even smaller delay time, of < 100 Myr (Man- 
nucci et al. 2005, 2006; Aubourg et al. 2008). The fact that 
our first time bin goes up to 420 Myr could then lead to a 
dilution of $1, due to its being averaged over the lower rates 
at these larger delays. 

Turning to the later bins in the recovered SN la DTD, 
the intermediate-delay bin of 420 Myr to 2.4 Gyr has a signal 
of * 2 = 0.0037tHoi4 SNe yr'^lO 10 Me)" 1 . Although the 
uncertainty is still relatively large, there is a ~ 2.5a DTD 
signal above zero in this bin. Naturally, some or all of this 
signal could be due to the "leak" from the prompt bin, dis- 
cussed above. Based on the observed 40% leak in the CC SN 
DTD, ~ 20% of the $2 signal could be leaked from $1. 

Finally, the last DTD bin, > 2.4 Gyr, has the clearest 
signal, #3 = (2.55t°;£j) x HT 4 SNe yr _1 (10 10 M ) _1 . At 
face value, this implies a ~ 4<r detection of a population 
of SNe la with large delays. The "A" parameter defined by 
Scannapieco & Bildsten (2005) was meant to measure the 
SN la rate per unit mass in an old population that has no 
ongoing star formation. Actually, SN la rates are expected 
to vary significantly among quiescent stellar populations of 
different ages, and therefore like the B parameter, A will 
depend on the ages of the stellar populations probed by a 
SN survey. Indeed, the typical values found for the A pa- 
rameter have been in the range A ~ (2 — 10) x 10~ 4 SNe 
yr _1 (10 10 M ) _1 (see compilation in Maoz 2008). However, 
^3 is the rate per unit mass formed. The published rates 
in old populations are a per unit of existing stellar mass 
(i.e., in stars and stellar remnants). The difference is given 
by the stellar mass returned to the interstellar medium via 
SN explosions and mass loss during stellar evolution. The 
fraction of returned material is an increasing function of 
time, and corresponds to ~25% after 0.1 Gyr and ~50% 
after 12 Gyr (Bruzual & Chariot 2003). Thus, for compar- 
ison to the A parameter, ^3 needs to be doubled, giving a 
rate of (5.1±H) x 10~ 4 SNe yr -1 (10 10 M ) -1 , in excellent 
agreement with other measurements of SN la rates in old 
populations. In principle, some or all of the ^ signal could 
again be the result of a leak from the prompt bin, due to the 
limited spatial coverage of the SDDS fibres. In practice, the 
time integral over the $3 rate, ^3At3, is seen to be compa- 
rable to ^lAti, as opposed to the few-percent leak into the 
* 3 bin that we found in the CC SN case in 33.4.21 This sug- 
gests that our ^3 value is real and largely uncontaminated. 
Furthermore, examining the time integrals over the best-fit 
DTD in each of the bins, and attempting to correct for the 
leak from bin 1 to bin 2, suggests a relative contribution 



to the total SN la numbers of (prompt:medium:delayed) sa 
2:2:1. However, this is subject to large statistical and sys- 
tematic uncertainties. 

It has been known for a long time that SNe la can occur 
in early-type galaxies with little or no star formation, and 
hence our measurement of a significant delayed component is 
hardly revolutionary. However, a SN la in a particular early- 
type galaxy can always be attributed to some residual low- 
level star formation combined with a short SN la delay time. 
Our measurement, on the other hand, provides a statistically 
robust determination of the delayed component, and its level 
relative to the prompt components, in a population with 
detailed measured SFHs. 

3.5.2 Full-sample SN la DTD, two-parameter dust model 

Using the second VESPA SFH reconstruction that utilises a 
two-parameter dust model, the best-fit DTD values (Fig. 2] 
and Table [1} for the prompt, medium, and delayed compo- 
nents are somewhat different from those obtained with the 
one-parameter dust model, However, these systematic dif- 
ferences are smaller than the statistical uncertainties of the 
results quoted above. For larger SN samples with smaller 
statistical uncertainties, the systematic uncertainties due to 
SFH modeling may, however, become dominant. In what fol- 
lows, we will use only the VESPA SFHs based on the first, 
one-parameter dust model. 

3.5.3 Culled-sample SN la DTD 

As in i)3.4.2l we have repeated the DTD derivation for SNe la 
from the subsample of 1900 galaxies that excludes Hubble 
types Sa to Sbc. This sample hosts 49 SNe la. Figure[5]shows 
the SN la DTD obtained for this culled sample. 

As was the case for the CC SNe in the culled sample, 
the prompt-bin amplitude is enhanced compared to the full 
sample DTD, and the medium-delay bin rate is lowered. The 
rate in the third bin, in contrast, is hardly changed compared 
to Figure [3] reconfirming the reality of that signal. 

3.5.4 Comparison to DTD models 

It is beyond the scope of this paper to attempt a detailed 
comparison of our recovered DTD to the many models that 
have been proposed, and to address the progenitor issue. 
Nonetheless, we have superimposed in Figure [5] which is our 
most reliable SN la DTD (insofar as the leakage problems 
are partly mitigated in it), a selection of DTDs that have 
appeared in the literature. Some are empirically motivated, 
in the sense that they were used to relate SN-rate data with 
star-formation measurements, using one of the methods out- 
lined in §1. Others are more theoretically motivated, based 
on a progenitor scenario. 

The Scannapieco & Bildsten (2005) "A + B" model has 
been discussed several times above. It is not exactly a DTD, 
but rather a prescription for relating a SN rate to a galaxy of 
a given current mass and SFR. Nevertheless, we have over- 
plotted this model in Figure [5] using the median literature 
values of the A and B parameters compiled by Maoz (2008) . 
We have associated the prompt B component with the same 
420 Myr bin of our DTD reconstructions, and the level to 
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Figure 5. Best-fit SN la delay-time distribution found for a sub- 
sample obtained by culling the LOSS-SDSS sample of all galaxies 
of Hubble type Sa-Sbc, for which SDSS fibre spectra do not rep- 
resent well the full stellar population of a galaxy. Symbols are as 
in Fig. [3] Note the enhanced ratio of rates between the first and 
second bins, due to reduced cross-talk. Also plotted, for compari- 
son, are several empirically or theoretically motivated DTD mod- 
els: A4-B, representation of the Scannapieco & Bildsten (2005) 
U A + B" model, plotted at the median levels of various esti- 
mates, as compiled by Maoz (2008); G, Gaussian DTD centred at 
4 Gyr with half-width 0.8 Gyr, proposed by Strolger et al. (2004); 
P, t~ - 5 power law, suggested by Pritchet et al. (2008); DD3, 
double-degenerate "DD-Close-3" model by Greggio (2005) with 
3Mq minimum initial mass; and SD, single-degenerate model of 
Greggio (2005), as shown by Greggio et al. (2008). The latter three 
models have been normalised to go through the best-fit observed 
LOSS-SDSS *3 rate for this subsample. 



be plotted is then just the median level, B — 2 x 10 3 Mg 1 
divided by this time interval. The median value for the A 
parameter is A = 6 x 1CT 4 SNe yr _1 (10 10 M Q ) _1 , which we 
plot as a constant rate at times > 420 Myr, but halved to 
account for stellar mass loss (see above). As already noted, 
"A + B" is an oversimplification, but its measured values are 
seen to be consistent with our DTD. 

Strolger et al. (2004) and Dahlen et al. (2004) have de- 
duced a Gaussian DTD centred at 4 Gyr with half-width 
0.8 Gyr, as a best fit to their comparison of SN rates to cos- 
mic SFH, out to z = 1.6. As seen in Figure [5] such a DTD 
is strongly at odds with our local, directly derived, DTD. 

Pritchet et al. (2008) have argued for a power-law DTD 
form, f- a - 5±0 - 2 ; based on analysis of the Supernova Legacy 
Survey data. A dependence of roughly t -0,5 is also expected 
simply from the formation rate of WDs, when considering 
the IMF and main-sequence stellar lifetimes. In Figure [5] 
we have plotted a t~ 0,5 curve, normalised to pass through 
the ^3 rate, which is the most robust and stable rate in our 
DTD reconstruction. It appears that a t~ ' 5 dependence is 
too shallow to match our derived DTD, especially consider- 
ing that 9\ is likely underestimated by us, given the possi- 
bility of remaining leakage from *t>i to ^2, even in our culled 
sample (as evidenced from the CC SN sample). 

Finally, we have overplotted in Figure [5] two of the an- 
alytical models of Greggio (2005), as presented by Greg- 
gio, Renzini, & Daddi(2008), based on stellar-evolution ar- 
guments and on various parametrisations of the possible re- 
sults of the complex common-envelope phases through which 



SN la progenitor systems must pass. For each of several 
SN la channels, Greggio (2005) calculates the DTDs that 
emerge when varying the values for the parameters describ- 
ing the initial conditions, and the mass and separation dis- 
tributions and limits of the systems that eventually explode. 
We show here one SD model and one DD model, again nor- 
malised to pass through our best-fit ^3 rate. 

The "DD-Close-3" label of the DD model refers to one of 
two possible parametric schemes used by Greggio (2005) to 
describe the WD separation distribution after the common- 
envelope phase, and to a minimum assumed initial mass of 
the secondary star in the binary, of 3Mq. This DD model 
appears to match well our recovered DTD, especially con- 
sidering that we have residual leakage problems from ^1 to 
^2, and hence *t>i is underestimated. It is easy to see from 
Figure [5] that this model, after its initial rise to maximum, 
is essentially a broken power law. At t < 400 Myr, the slope 
is —0.5, just like the Pritchet et al. (2008) model slope. This 
slope is the result of stellar-evolution lifetimes and IMF (see 
above). At t > 400 Myr, the slope is —1.3. A slope of roughly 
— 1 is generic to models in which the merger rate is deter- 
mined by energy loss to gravitational waves (e.g., Greggio 
2005; Totani et al. 2008). 

The Greggio (2005) SD model also matches our re- 
covered LOSS-SDSS DTD in Figure [5] remarkably well. It 
should be remembered, however, that in a full physical 
model, the normalisation is not a free parameter, and is dic- 
tated by the efficiency of SN la production from the poten- 
tial progenitor population. As emphasised by Maoz (2008), 
and further discussed in £13.6.21 below, the efficiencies found 
by binary population synthesis models are at least a fac- 
tor of a few, and likely an order of magnitude, lower than 
indicated by observed SN rates (see also Mennekens et al. 
2010). Furthermore, the SD model, whose shape matches 
our derived DTD so nicely, has been found to be even less 
efficient, by yet another order of magnitude, in some studies 
(e.g., Tutukov & Yungelson 2002). 

From this brief comparison of our recovered DTD to 
previous work, we conclude the following. (1) The monoton- 
ically decreasing nature of the DTD that we have found is 
in agreement with most previous empirical determinations 
and theoretical expectations, except for that of Strolger et 
al. (2004). (2) The decline with time of our DTD is steeper 
than in the Pritchet et al. (2008) model. Two of the Greggio 
(2005) models, SD and DD-Close-3, fit well the shape of the 
recovered LOSS-SDSS DTD. (3) The normalisation of full 
physical models to the observations has not been considered 
here, but it is another point at which models are challenged 
by these and other observations of SN rates. 



3.6 The time-integrated DTD 

Another interesting observable is the integral of the DTD 
over the age of the universe (to), 



N SN /M = #(t) dt 



(9) 



which gives the total number of SNe that eventually explode, 
per unit stellar mass formed in a short burst of star forma- 
tion. We can approximate the integral in Eq. [9] with a sum 
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over the binned DTD, 



i=l,K 



(10) 



The adequacy of this approximation will depend on the true 
form of i&(t) and on the number of bins. From the simula- 
tions described in the next section, we find that, for de- 
clining power-law DTDs represented with three time bins, 
Eq. [TU] typically underestimates JVsn /M systematically by 
~ 10 - 20%. 



3.6.1 Core-collapse SN yield per stellar mass 

If all stars with mass > 8 Mq explode as CC SNe, then 
for CC SNe, Nsn/M is just the ratio of the number of 
stars formed above this mass limit to the total stellar 
mass formed, and is easily calculated to be Ncc/M = 
0.010 SNe Mq 1 for the "diet" Salpeter IMF (as well as for the 
Kroupa (2007) IMF assumed by the the VESPA SFH recon- 
struction). In contrast, the integral over our reconstructed 
DTD of the CC SNe is (0.00585 ± 0.00085) SNe Mq 1 , a fac- 
tor of ~ 2 lower than expectecjf|. However, we can show that 
this deficit of CC SNe largely disappears in various subsam- 
ples of the LOSS-SDSS sample, and is therefore not a real 
effect. 

For example, we have reduced the LOSS-SDSS galaxy 
sample by limiting the maximum galaxy distance to progres- 
sively smaller values. For each distance-limited sample, we 
recover the DTD. Figure [5] shows the values of Ncc/M we 
obtain by integrating over the reconstructed CC SN DTDs, 
for each distance limit. Several of these values are also listed 
in Table [T] Clearly, Ncc/M rises with decreasing distance, 
and reaches close to the theoretically expected value when 
the sample is limited to galaxies closer than D ~ 65 Mpc. 
From this plot we estimate the observed limiting value to be 
Ncc/M = 0.010 ±0.002. 

Although a potential explanation for the dependence of 
Ncc /M on sample volume would be an increasing fraction of 
missing CC SNe with distance, perhaps due to mild obscu- 
ration of such SNe by dust, this is unlikely to be the correct 
explanation. Ncc/M is derived from the DTD, which takes 
into account the visibility time of each type of CC SN at the 
distance of each galaxy in the sample. The visibility time, in 
turn, was calculated, as described in Li et al. (2010a), using 
the observed LF of 87 nearby CC SNe, in which the effects 
of extinction are automatically included. 

More likely, the decrease of Ncc/M with distance is 
another manifestation of the SDSS fibre aperture problem. 
At larger distances, there are fewer low-mass galaxies and 
late-type spirals in the LOSS sample, and more high-mass 
intermediate and early types. In the earlier-type galaxies, 
VESPA in more likely to overestimate the total galaxy mass, 
based on the spectrum of the bulge population probed by 
the SDSS aperture. Indeed, if we rederive the CC SN DTD 
and Ncc/M, but exclude the earlier Hubble types (as in 
33.4,2( 1. or exclude massive galaxies, the values of Ncc/M 
for the samples limited to D < 125 Mpc, D < 150 Mpc, 



3 Here, the statistical error on Ncc/M is found using the covari- 
ance matrix [C] for "J/, and the vector of time-bin intervals At = 
(Ati, At 2 ,...,At K ), such that A(JV S n/M) = (At[C] At) 1 ^. 
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Figure 6. The time-integrated DTD, i.e., the total number of 
SNe produced over a Hubble time, per unit stellar mass formed, 
for different distance-limited galaxy subsamples. Top full symbols 
are for CC SNe (SNe II plus SNe Ibc), bottom full symbols are 
for SNe la. Top empty symbols are for CC SNe in a subsample of 
galaxies with masses < 3 X 10 10 Mq. Bottom empty symbols are 
for SNe la in a subsample of galaxies with masses < 7 X 10 10 Mq . 
The rise in Ncc/M toward its theoretically expected value with 
decreasing maximum sample distance or with sample culling is 
another manifestation of the SDSS small-fibre-aperture problem. 
The value of Ncc/M at < 65 Mpc indicates that most stars 
above 8 Mq produce CC SNe. The constant value of 7Vi a /M, on 
the other hand, shows that the systematics affecting SNe la are 
less dependent on distance. Nevertheless, culling the higher-mass 
galaxies from the sample increases N\ & /M at all distances. The 
time-integrated ratio of CC SNe to SNe la is seen to be roughly 
4:1. 



and D < 200 Mpc rise significantly, and approach their 
D < 100 Mpc values. For example, we show in Figure [6] 
(upper empty symbols) that, if we limit the sample to galax- 
ies less massive (based on their VESPA reconstructions) 
than 3 x 10 10 Mq, the best-fit DTDs at those distances give 
Ncc/M = 0.0086 ± 0.0017, much closer to the expected 
value. (Naturally, the reduced sample sizes lead to larger 
statistical error bars in Fig. [6]) Another indication of the 
nature of the effect can be seen in Table 1, where the ra- 
tio of $2/^1) which quantifies the leak from the first to the 
second bin, decreases for the CC SNe as D is reduced. 

The fact that a galaxy sample can be defined (the 
D < 65 Mpc sample in the above example), for which 
Ncc/M actually reaches its expected value, is important. 
It confirms that, indeed, the majority of stars with > 8Mq 
produce CC SNe. This cannot be taken for granted. Given 
current observational and theoretical limits, the low-mass 
limit for core collapse could be as high as 10 or even 11 Mq 
(see Smartt et al. 2009, and references therein). Moving 
the limit from 8 to 10 Mq would decrease the expected 
Ncc/M by about 30% (i.e., to our 2a observed lower limit 
on Ncc/M). Raskin et al. (2008) have recently obtained a 
similar estimate of the lower mass limit for core collapse 
by matching, on the one hand, the differences in the spa- 
tial distributions of stars and SNe in spiral galaxies to, on 
the other hand, the predictions of simple stellar-population 
aging models. 

Furthermore, in principle, a sizable fraction of high- 
mass stars could collapse directly to a black hole, with 
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only a weak or null SN explosion accompanying the col- 
lapse (Kochanek et al. 2008). Our measurements of Ncc/M 
also argue against this scenario, unless the minimum mass 
for CC SNe is even lower than 8 Mg . We note that these 
conclusions are reinforced by the fact that the discrete sum 
in Eq. 1101 underestimates the true integral over 

3.6.2 The SN la yield per stellar mass 

Figure [6] also shows the results of a similar analysis for the 
SNe la: the integral Ni & /M over the best-fit SN la DTD 
(shown in Fig. [3] for the full sample), again as a function of 
the sample distance limit. We see (lower filled symbols) that 
there is no obvious dependence on distance, with Ni^/M ap- 
parently constant at ~ 0.00140 ± 0.00033 SNe Mg 1 . How- 
ever, if we again limit the sample in mass, for example to 
< 7 x 1O 1O M0, we obtain (lower empty symbols in Fig. 
higher values, iVi a /M w 0.00230±0.00060 Mg 1 . The system- 
atics of the masses of the SN la host galaxies, while present, 
are apparently less dependent on distance than those of the 
CC SN hosts. 

These values of Ni^/M are also of interest. As al- 
ready noted above, previous studies have found, just for 
the prompt SN la component embodied in the B param- 
eter, values of B = (1 - 3) X 10~ 3 Mq\ We note that the B 
parameter, like Ni^/M, relates to the stellar mass formed. 
Manmicci et al. (2006) obtained an empirical DTD based on 
the observations described by Mannucci et al. (2005) and 
Delia Valle et al. (2005). Their time-integrated SN rate is 
0.0013 SNe Mg 1 . However, this rate is per unit of existing 
stellar mass, after assuming a fraction of recycled gas of 0.30 
as an average among populations of different ages (Bruzual 
& Chariot 2003). Accounting for this factor, the Mannucci 
et al. (2006) result corresponds to ~ 0.0010 SNe Mg 1 of 
formed stars, similar to the values measured here. 

As discussed at length by Maoz (2008), Ai a /M is one 
of several observables that can be related directly to the 
fraction r\ of stars in some initial mass range [mi, 7712] that 
eventually explode as SNe la: 

_ Ma JoT m{dN / dm)dm 
77 ~ M I^(dN/dm)dm ' 1 ' 

where dN/dm is the IMF. For the "diet Salpeter" IMF 
(which, again, gives results similar to the Kroupa IMF as- 
sumed by the VESPA SFH reconstruction), and an initial 
mass range of 3-8 Mg, often considered for the primary 
stars of SN la progenitor systems, the ratio of the two in- 
tegrals equals 33. Adopting the higher value of Ni^/M that 
we have found in the culled sample, under the assumption 
that it is more robust against the fibre-aperture effect, we 
obtain 77 = 7.6 ± 2.0%. This is in agreement with the re- 
sults of Mannucci et al. (2006), who found 77 = 4.3% for 
these parameters. As noted by Maoz (2008), the consistently 
high values of the exploding fraction, 77, derived from several 
different observables, may constitute a problem for current 
progenitor models. 

Interestingly, in galaxy clusters, Maoz et al. (2010) es- 
timate a time-integrated number of SNe la per present-day 
stellar mass of 0.011 SNe Mg 1 . This estimate is based on 
the observed ratio of total mass of iron (both in stars and 
in the intracluster medium) and the mass in stars. For an 



IMF with a standard high-mass end, the present-day mass 
observed in stars is related to the number of CC SNe, whose 
contribution to the iron mass has been estimated and sub- 
tracted. The remaining iron mass is then due to the SNe la. 
Multiplying by 0.5 to convert to formed mass, rather than 
present mass, this gives Ni^/M = 0.006 SNe Mg . Maoz 
et al. (2010) show that, given the uncertainties in observed 
cluster properties, this number could decrease by perhaps a 
factor of 0.6 at most, to Nj^/M = 0.0035 SNe Mg 1 . This 
is still at least a factor of 1.5 greater than, but marginally 
consistent with, the N u /M = 0.0023±0.0006 SNe Mg 1 that 
we have found here based on LOSS. 

The cluster-based value of A r i a /M may be evidence for 
early enrichment of clusters by CC SNe from a top-heavy 
IMF. These CC SNe from massive stars, which left no traces 
in the form of low-mass relatives, would have then produced 
the bulk of the iron mass in clusters. Alternatively, the large 
iron mass in clusters could have indeed come from SNe la, 
but this would imply a more efficient production of SNe la 
in cluster environments. Intriguingly, Sharon et al. (2007), 
Mannucci et al. (2008), and Graham et al. (2008) have all 
found evidence for SN la rates enhanced by factors of a few 
in cluster galaxies, compared to field early-type galaxies. 



3.6.3 The ratio of CC SNe to SNe la 

From the ratio of the time-integrated DTDs, Ncc/M and 
Nu/M (using the value of Ncc/M at D < 70 Mpc for 
which the CC SN counts are fairly complete, but for which 
the errors are not excessively large because of the limited 
sample size, and Ni^/M from the low-mass sample, see 
above), the time-integrated ratio of CC SNe to SNe la is 
Ncc/Nu, = 4l 3 .5- If we force the value of Ncc/M = 0.01 
expected from a diet-Salpeter IMF with a low-mass CC 
limit of 8 Mg, the allowed range in the ratio shrinks to 
Ncc/Nia, = 4^q 7 ). We note that the ratio of time-integrated 
DTDs is distinct from the observed ratio of current rates, 
which is measured to be about 3:1 in local surveys (Man- 
nucci et al. 2005; Li et al. 2010a). The ratio of current rates 
depends on the summed SFH of the galaxies in the volume. 
It can be arbitrarily high for very young, star-forming pop- 
ulations (in which few SNe la have had time to form yet) 
to zero for old, inactive populations (with no CC SNe). In 
contrast, the ratio of the time-integrated DTDs, like the 
DTDs themselves, is independent of SFH, and intrinsic to 
the stellar-evolution processes that lead to SNe. 

The ratio Ace /Ma of course equals 
(0.01Mg 1 )/(Ai a /M), for the diet Salpeter IMF and 
the said lower limit for CC SNe, and thus the cluster-based 
lower limit of Ai a /M > 0.0035 SNe M" 1 (Maoz et al. 2010), 
discussed above, implies Acc/Ai a < 2.9, and possibly 
even 1:1. De Plaa et al. (2007), comparing cluster element 
abundances to theoretical SN element yields, have deduced 
an integrated CC SN to SN la ratio of 1:1. On the other 
hand, the observational uncertainties are such that a time- 
integrated CC SN to SN la ratio of roughly 3:1, despite 
some tension, is consistent both with cluster measurements 
and with the LOSS data we have analysed here. 
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4 TESTS ON SIMULATED SAMPLES 

To test our DTD recovery procedure and examine its per- 
formance on different kinds of input datasets, we have re- 
peated the Monte-Carlo mock survey generation described 
in §2, above, using the LOSS-SDSS galaxy sample and its 
VESPA-based SFHs (with one dust parameter). However, 
rather than inputing the best-fit DTD from the inversion of 
the real data, we can choose other model DTDs and convolve 
those DTDs with rxiij according to Eq.[2]-[3] or we can change 
the sample properties in additional ways, as described below. 
Using the actual LOSS galaxy SFHs and visibility times (as 
opposed to, for example, random SFHs and visibility times) 
makes for a more realistic simulation. As before, the number 
of SNe found in each galaxy in each mock survey is drawn 
from a Poisson distribution with expectation value Ai. The 
different-shaped DTDs which can be input and recovered 
can also be scaled up or down to produce a larger or smaller 
total number of SNe in the mock survey (or, equivalently, 
the visibility times, U, can be scaled up or down). In these 
simulations, the SFHs are assumed to be error-free inde- 
pendent variables, and are used both for creating the mock 
samples and deriving their DTD. Therefore, these simula- 
tions will not display the "leakage" between bins that we 
have encountered with the LOSS-SDSS sample, nor other 
problems due to systematic or random errors in the SFHs 
of the galaxies. 

From our simulations, we find that the distribution of 
output DTD amplitudes is centered on the input DTD val- 
ues, meaning that the method reliably recovers the input 
DTD with little bias. Nonetheless, for low DTD amplitudes 
combined with large uncertainties (due to small SN num- 
bers, see below), there can be some "pile-up" in the distri- 
bution at zero amplitude, as a result of the positivity con- 
straint of the DTD. Thus, obtaining a zero amplitude for a 
bin in the DTD reconstruction can happen, even when in 
reality the amplitude is nonzero but low, and the number of 
SNe in the survey is small. 

We find that the relative uncertainty in the DTD am- 
plitude in the jth bin scales roughly as one would expect 
from Poisson statistics, considering the number of SNe that 
contribute to every bin in the DTD. For example, with three 
time bins, the relative error in the first bin of the DTD is 



A/tot 



Ei + Ei m.,2*2 + Ei "»i,3*; 



-1/2 



(12) 

For a survey with a fixed total number of SNe, N to t, there 
will thus be a tradeoff in the analysis between DTD accu- 
racy and resolution. The uncertainties do not depend on 
the number of galaxies monitored, and thus a brief-duration 
survey of many galaxies and a long-duration survey of few 
galaxies are equivalent, as long as they produce the same 
total number of SNe. 

Figure [7] examines the quality of the reconstruction as 
the number of time bins is varied between three and five. 
Here, we have assumed a survey with 15,000 galaxies and 
~ 370 SNe, similar to the full LOSS SN la sample, if all 
galaxies in it had SFH reconstructions available. To pro- 
duce this large mock sample, we simply clone several times 
the 3505 SFHs and visibility times of the LOSS-SDSS sam- 
ple, with its three time bins. To obtain a sample with five 
time bins in each galaxy's SFH, we have split, in this exam- 
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Figure 7. DTD recovered using three or five bins from a mock 
sample of 15,000 galaxies hosting 370 SNe la, the numbers char- 
acterising the full LOSS sample. The input DTD is a t^ 1 power 
law (dashed line). Points are plotted at the value of the input 
binned DTD, and error bars show the ranges that include 68% of 
the recovered best-fit solutions in repeated random realizations 
of SN surveys on the mock sample. 



pie, the first SFH bin into two bins: a 0-100 Myr bin and a 
100-420 Myr bin. The stellar mass in each sub-bin was ran- 
domised by ±50% around the mass value corresponding to 
the relative time fractions. (Such randomisation is essential, 
as otherwise the rriij are no longer independent variables.). 
Similarly, we further split the last bin into two bins, corre- 
sponding to 2.4-6.5 Gyr and 6.5-13.5 Gyr. For the DTD in 
these simulations, we have taken a ^(t) oc t -1 dependence 
(which is roughly generic to DD models; see t|3.5.4l above), 
scaled so as to give the desired total number of SNe using 
the LOSS visibility times. The simulations show that surveys 
with several hundred SNe enable reasonably good DTD re- 
covery in terms of both accuracy and temporal resolution. 
Specifically, Figure [7] demonstrates that already-existing SN 
surveys, and LOSS in particular, have the power to measure 
reliably the SN la DTD, if the SFHs could be estimated in 
a comprehensive and unbiased way for the full sample of 
~ 15, 000 galaxies. We briefly discuss the prospects for this 
in the concluding section, below. 

We have gauged the effect of the time binning on the 
time integral over the DTD, Nsn/M, as approximated in 
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Eg. 1101 To do this, we create mock samples using large num- 
bers of bins in the SFH and the DTD. The mock samples also 
have a large number of SNe in order to minimise the statisti- 
cal error and isolate the systematic effect due to the binning. 
We then rebin the SFHs of the sample into the same three 
coarse time bins we have used for the LOSS-SDSS analy- 
sis, and we recover the DTD binned into those same three 
time intervals. Finally, we compare the integrals Nsn/M 
over the input DTD and the recovered DTD. We find that 
the iVsN /M of the recovered DTD systematically underesti- 
mates the input Nsn/M by 10-20%, with the larger values 
for input DTDs that rise more steeply at small delays. This 
systematic error is comparable to the statistical errors in 
Nsn/M we have found for the LOSS-DTD sample. Larger 
SN samples will naturally allow finer temporal binning of 
the DTD, reducing this systematic effect. 



5 CONCLUSIONS AND OUTLOOK 

We have presented a method to recover the SN DTD by 
jointly analysing SN survey data and reconstructed SFHs 
of the individual galaxies in a survey, while accounting for 
the small-number, Poisson-statistical nature of SN events in 
those individual galaxies. Our method is an improvement 
over previous ones in that it uses the full information con- 
tained in the data and avoids unnecessary averaging. The 
method is based on casting the expected SN numbers in 
each galaxy in the survey as a set of linear equations, in 
which the parameters to be constrained by the data are the 
values of the DTD in binned time intervals. We have shown a 
maximum-likelihood method by which to invert those equa- 
tions to recover the DTD, and demonstrated the method's 
performance using simulated mock samples. We have ap- 
plied the method to a sample consisting of those galaxies 
in the LOSS SN survey that have SDSS-VESPA SFH re- 
constructions, along with the SNe that these galaxies have 
hosted during the LOSS. 

The recovered DTDs for CC SNe and SNe la are limited 
by Poisson statistics, because of the small numbers of SNe 
in this restricted subset of LOSS, and by systematic errors, 
because for many galaxies the light within the SDSS fibre 
aperture does not reliably represent the full stellar popula- 
tion of the galaxy and its SFH. Despite these shortcomings, 
we are able to make the following statements based on anal- 
ysis of these data. 

(1) A "prompt" SN la population, defined here as one that 
explodes within 420 Myr of star formation, exists at > 99% 
confidence. This confirms the growing number of reports of 
such a population, at levels similar to those found here. 

(2) We clearly detect, at the 4<j level, a "delayed" SN la 
population with delays in the range 2.4-13 Gyr, and a mean 
rate over this interval of * 3 = (2.6 ± 0.7) x 10~ 4 SNe yr" 1 
per 1O 1O M formed, or A = (5.1 ± 1.4) x 10~ 4 SNe yr" 1 
per remaining 10 10 M in an old population. These first two 
conclusions together indicate that the SN la DTD peaks at 
short delays, but extends over a broad range of delay times, 
out to at least several Gyr. Progenitor models, a few of which 
we have briefly compared to our recovered DTD, will need 
to reproduce the observed numbers. 

(3) The time-integrated SN la yield is Nu/M = (2.3 ±0.6) x 
10~ 3 SNe per unit solar mass formed, or (4.6 ± 1.2) x 10~ 3 



SNe la per remaining unit solar mass in an old population. 
The best-fit value is a factor of 1.5-3 lower than the corre- 
sponding number in galaxy clusters, as deduced from their 
measured iron-to-stellar mass ratios. Although the current 
uncertainties can still accommodate this difference, it may 
indicate that clusters underwent additional enrichment by 
CC SNe from an early stellar population with a top-heavy 
IMF, or that SN la production is more efficient in galaxy 
clusters than in the field. The latter possibility has support 
both from direct cluster SN rate measurements and from 
cluster element abundance analysis. The measured time- 
integrated SN la yield also implies r\ ~ 8% for the exploding 
fraction among the parent population of SN la primaries, if 
assumed to come from the 3-8 Mq initial mass range. 

(4) The time-integrated CC SN yield is N C c/M = (1.0 ± 
0.2) x 10~ 2 SNe per unit solar mass formed. This rules 
out low-mass limits for CC explosions that are much above 
8 Mq . Conversely, scenarios in which a significant fraction of 
high-mass stars end their evolution without SN explosions 
are excluded, unless the low-mass limit for core collapse is 
significantly below 8 Mq. 

(5) The ratio of CC SN to SN la numbers from a brief burst 
of star formation, integrated over a Hubble time, is (4 jl i 5). 

Our work points to the kinds of data that would improve 
upon these results. First, a larger sample of survey galaxies 
with spectroscopy, and hence SFHs, would obviously reduce 
the Poisson errors and permit better temporal resolution. 
As there are 14882 galaxies monitored by LOSS, obtaining 
spectra for most or all of them would be a large, but not 
impossible, task. A dedicated or partly dedicated 2-4 m- 
class telescope could achieve this on a few-year timescale. 
Long-slit spectra, drifted across the galaxy perpendicular to 
the slit length, would be more representative of the entire 
stellar population of each galaxy, largely avoiding the small- 
fibre-aperture problem we have encountered. Ideally, instead 
of long-slit spectra of the LOSS galaxies, one would obtain 
integral-field spectroscopy of each of these galaxies (using, 
e.g., an instrument such as SAURON on the William Her- 
schel 4.2 m telescope; Bacon et al. 2001). In addition to 
including the full SFH of each galaxy without any aperture 
losses, such data, in the context of our method, would easily 
allow breaking up each galaxy into independent subunits, 
and considering the SFH of each subunit in relation to the 
SNe that it hosted (or, for almost all subunits, the SNe that 
it did not host). 

It might be objected at this point, that because of the 
random stellar velocities in each galaxy, SN progenitor stars 
diffuse away from the regions where they were formed, and 
that this would invalidate our approach; by the time a SN la 
exploded, it would reside within a region that is completely 
different from the one in which its progenitor was formed. 
While this is true, and it would in fact invalidate a tradi- 
tional "SN delay time" analysis in which a single charac- 
teristic stellar population age is assigned to a region, it is 
inconsequential to our current method, in which the entire 
SFH of the region is considered. The reason is that the same 
spatial diffusion affects both the progenitor population and 
those stars among it that eventually explode. To see this, 
consider, as a toy example, a grid of 3 x 3 adjacent "cells" 
in a particular galaxy. Suppose, for example, that 500 Myr 
ago there was a short burst of star formation in the central 
cell, forming a stellar mass M, and no activity in the other 
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cells. Suppose, further, that the SN DTD is such that the 
stellar population formed in the burst leads, 500 Myr later, 
to nine SNe over the course of a decade, which are there- 
fore detected in this galaxy by a SN survey such as LOSS 
(this is admittedly a somewhat unrealistically large number 
of SNe, but is used here just for the sake of illustration). 
The galaxy has thus produced a ratio of 9/M SNe per unit 
stellar mass formed 500 Myr ago. Finally, suppose that the 
stellar diffusion timescale in the galaxy is such that, over 
the 500 Myr, the progenitors of the 9 SNe, before explod- 
ing, have drifted out of the central cell in which they were 
formed, and there is now, on average, one SN in each cell. 
However, the entire stellar population of the burst will have 
diffused in the same way, and therefore each cell will have 
1/9 of the 500-Myr-old population that was originally in the 
central cell. When we compare SN numbers to the 500 Myr- 
old stellar mass present today in each cell, we will see 1 SN 
per M/9 of stellar mass formed. In the DTD we derive, we 
will therefore still deduce the correct ratio of 9/M SNe per 
unit stellar mass formed 500 Myr earlier. 

This argument holds, no matter what are the lookback 
times or the diffusion timescales. It also holds for arbitrarily 
complex SFHs, which can be viewed as linear superpositions 
in space and in time of toy models of the type above. If a 
past starburst at a certain place and time in a galaxy pro- 
duces SNe that are observed in the course of a survey, the 
stellar population of the burst and the SNe that it produces 
will both diffuse in the same way. If, on the other hand, a 
specific burst does not produce SNe detected by the survey 
(because the DTD has a low amplitude at the corresponding 
delay), then there will be no correlation between the number 
of SNe per cell and the mass of stars of that age per cell. 
An individual cell hosting SNe may, of course, include un- 
related stellar populations that did not produce those SNe, 
but whose stars nonetheless drifted into the cell. However, 
over the entire galaxy, there will be no correlation between 
SNe and stars of that particular age, and it is such correla- 
tions that drive the results of our DTD recovery method. 

We note that for the integral-field spectroscopy ap- 
proach to work, the signal-to-noise ratio of the spectra of 
each individual galaxy cell needs to be sufficiently high for 
a reliable SFH reconstruction to be obtained. In partic- 
ular, the presence of old stellar populations that are su- 
perimposed on younger and more luminous stars must be 
detectable. Naturally, spatially resolved, medium-spectral- 
resolution data of such a large sample of nearby galaxies 
would find many additional applications, and hence such 
data are worth the large effort required. 

Shortly before submission of this work, Brandt et al. 
(2010) presented a DTD reconstruction analysis of a differ- 
ent SN sample. Their methodology shares several elements 
with ours. Brandt et al. (2010) study 107 SNe la from SDSS- 
II. Like us, they use VESPA to derive SFHs for a sample of 
SDSS galaxies, binned into three time bins, identical to those 
we have chosen. Like us, they treat the DTD amplitudes in 
the three discrete bins as free parameters, which are deter- 
mined by a maximum-likelihood procedure. However, rather 
than comparing directly the presence or absence of SNe in 
each galaxy to the predictions of the DTD model (as we 
have), they use the DTD to create mock SN-host samples, 
and compare the mean spectrum of the mock host samples 
to the mean spectrum of the real host galaxies. Brandt et 



al. (2010) reach similar conclusions to ours, namely, signif- 
icant detections of both prompt (< 420 Myr) and delayed 
(> 2.4 Gyr) SN la DTD components. 
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